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Abstract 

We show that a non-equilibrium diffusive dynamics in a finite-dimensional space takes in the Lagrangian 
frame of its mean local velocity an equilibrium form with the detailed balance property. This explains the 
equilibrium nature of the fluctuation-dissipation relations in that frame observed previously. The general 
considerations are illustrated on few examples of stochastic particle dynamics. 



1 Introduction 

In the last decades, non-equilibrium statistical mechanics has been a subject of intensive studies. One 
of the multiple aims of the research is the understanding of essential differences between the equilibrium 
and non-equilibrium dynamics. This is the question that we shall address below. In the modelling 
of statistical-mechanical dynamics, an important role has been played by stochastic Markov processes. 
Although largely idealized, they often provide a sufficiently realistic description of experimental situations 
and have traditionally served as a playground for both theoretical considerations and numerical studies. 
The Markov processes corresponding to the equilibrium dynamics are characterized by the detailed balance 
property assuring that the net probability fluxes between micro-states of the system vanish. On the other 
hand, in the non-equilibrium Markov dynamics, the detailed balance is broken and there are non-zero 
probability fluxes even in a stationary situation. 

In the present paper, we shall consider only diffusive processes, discarding Markov processes with 
discrete time or random jumps. For such systems, the detailed balance can be expressed as the vanishing 
of the probability current that is non-zero in the non-equilibrium situations. It is convenient to represent 
the probability current in a hydrodynamical form as the instantaneous probability density of the process 
multiplied by the mean local velocity. The latter is the average instantaneous velocity of the process 
conditioned to pass through a given point. It will play the main role in what follows. 

In the past, there have been many attempts to apply ideas from statistical mechanics to the hydrody- 
namics of turbulent flows. The success was limited by the fact that most methods of statistical physics 
had been developed for systems in or close to equilibrium whereas developed turbulence is a far-from- 
equilibrium phenomenon. Here we shall follow a reversed strategy, applying an idea from hydrodynamics 
to non- equilibrium statistical mechanics. There is a long tradition (going back to Lagrange) to describe 
the evolution of hydrodynamical fields in the Lagrangian frame that moves with fiuid particles ^T\. It 
is believed that such a description makes the intrinsic features of fiuid dynamics at small scales more 
directly accessible than in the Eulerian (i.e. laboratory) frame. This is particularly true about the hy- 
drodynamical advection that gains a simple representation in the Lagrangian frame. The main result 
of the present paper consists of a simple observation that the non-vanishing probability 
current in a Markov diffusion may be decoupled from the stochastic dynamics by passing to 
the Lagrangian frame of the mean local velocity. More exactly, in the latter frame, the stochastic 
dynamics, although non-stationary, satisfies the detailed balance condition and the instantaneous proba- 
bility density of the process does not change in time. The equilibrium-like Lagrangian-frame process does 
not contain information about the non-vanishing probability current of the original Eulerian-frame pro- 
cess but, if that information is provided independently, the Eulerian-frame process may be reconstructed 
from the Lagrangian-frame one. In short, the passage to the Lagrangian frame of the mean local velocity 
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re-expresses a non-equilibrium diffusion process as an equilibrium-type one plus the decoupled probabil- 
ity current. To our knowledge, this rather straightforward observation about non-equilibrium diffusions 
has not been discussed in the literature, although a similar idea was recently employed in the quantum 
many-body dynamics [32] . 

The paper consists of seven Sections and four Appendices. Sect .[5] sets the stage and notations by 
briefly stating the basic definitions relevant for the diffusion processes that we consider. We introduce 
the notions of the probability current and of the mean local velocity and recall the concept of detailed 
balance. The crucial Sect. |3] is devoted to the Lagrangian picture of diffusions. We define the Lagrangian 
frame of the mean local velocity and compute the instantaneous probability density of the Lagrangian- 
frame process. By working out the stochastic differential equation satisfied by this process, we show 
that it is a non-stationary diffusion with the detailed balance property. Two simple examples illustrate 
the general considerations: a diffusion of a particle on a circle in the presence of a constant force and 
a linear stochastic equation describing a Rouse model of a polymer in shear flows. We also discuss the 
reconstruction of the original Eulerian-frame process from the Lagrangian-frame one. Sect.|3]is devoted to 
the Langevin equations with both Hamiltonian and non-conservative forces. In this case, it is convenient 
to modify the definition of the probability current and the mean local velocity to assure that they vanish 
in the absence of the non-conservative drift. The main properties of the Lagrangian-frame process are 
unaffected by this modification. We illustrate the general discussion by the example of a harmonic chain. 
Sect. l5.2l discusses the extensions of the Fluctuation-Dissipation Theorem to the non-equilibrium situation 
in the light of the results about the Lagrangian-frame process. These results provide a deeper reason for 
the observation made in [6], see also [28) . that the fiuctuation-dissipation relations takes the equilibrium 
form in the Lagrangian frame of the mean local velocity. In Sect.[6l we point out that important non- 
equilibrium diffusion processes in infinite-dimensional spaces, like the one-dimensional KPZ equation or 
the processes describing the large-deviations regime of fiuctuations around the hydrodynamical limit of 
the boundary-driven zero-range particle processes do not possess Lagrangian picture. Finally, Sect.[7] 
presents our conclusions. Appendices collect some more technical arguments. 

Acknowledgements. The authors thank Gregory Falkovich for discussions. R.C. acknowledge the 
support of the Koshland Center for Basic Research and K.G. of the project ANR-05-BLAN-0029-03. 



2 Eulerian picture of diffusions 

2.1 Diffusion processes 

We shall begin by considering a general diffusion process xt in a d-dimensional (phase-)space X with 
coordinates (a;'), of the same type as in ref. [5] that was devoted to the study of fiuctuation relations for 
such processes. The examples we shall have in mind include various types of Langevin dynamics used to 
model equilibrium and non-equilibrium dynamics as well as the Kraichnan model of turbulent advection 
|10) . Of the rich theory of diffusion processes, see e.g. [2611241131] . we shall need only few basic facts that 
we collect below. The process Xt is assumed to satisfy the stochastic differential equation (SDE) 

Xt = ut{xt) + Ct{xt) , (2.1) 

where it = and, on the right hand side, ut{x) is a time-dependent deterministic vector field (the 
drift), and Cti^) is a Gaussian random vector field with mean zero and covariance 

{Cdx)Ciy)) = 2Sit-s)Dnx,y). (2.2) 

Note that (t{x) is a white noise in time so that Eq. (|2.2|l requires a choice of a stochastic convention. As 
in [5], we shall interpret it in the Stratonovich sense to assure that ul{x) and Q{x) transform as vector 
fields under a change of coordinatei|3- The single time expectations of functions of the process Xt evolve 
according to the equation 

iifi^t)) ^ (iLtf){xt)), where Lt = uld, + d^d^' (2.3) 
at 

with 

d-Yix) = Dl^{x,x) ul{x) = ul{x) - rl{x) , rl{x) = dyj D'/ {x,y)\y=x (2.4) 
^In probabilists' notations. Eg. 112.11 1 would read dxt = ut(xt) dt + Y^, X„ (xt ) o dW^^ where Xn are vector fields such 

n 

that 2D^^ {x,y) = X^{x)X^{y) and W" are independent Wiener processes. 
n 
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are the instantaneous generators of the process xt- Note the presence of the term rt correcting the 
drift and due to the dependence the covariance of C,t on the points in X. The time evolution of the 
instantaneous (i.e. single-time) probability density function (PDF) of the process 

pt{x) = {6{x-xt)) (2.5) 

is governed by the formal adjoints lJ of the generators Lt : 

dtpt = Llpt = -di [u'tpt ~ d^t djpt] ■ (2.6) 

The transition PDF's of the Markov process xt given by the conditional expectations 

P{s,x;t,y) = ps(xy'^ {5{x ~ Xs)5{y ~ Xt)) (2.7) 

with s < t satisfy the Chapman-Kolmogorov composition rule J P{r,x; s, y) P(s, y; t, z) dy = P{r, x; t, z) 
and the Kolmogorov differential equations 

dsP{s,x-t,y) = -Ls{x)P{s,x-t,y), dtP{s,x-t,y) = Ll{y) P{s,x-t,y) . (2.8) 

The latter, together with the condition P{t,x;t,y) — S{x — y), determine the transition probabilities 
under appropriate regularity assumptions [31j . 

2.2 Probability current and mean local velocity 

Some other basic notions concerning Markov diffusions will play a central role below. The evolution 
equation (|2.6|l for the instantaneous PDF (|2.5p of the process Xt has the form of the continuity equation 

dtPt + V-jt = (2.9) 

with the probability current 

jl = [ul - d^d,]pt (2.10) 

whose flux through the boundary of any region V gives the rate of change of the probability that xt 
belongs to V. A more transparent interpretation of the current jt{x) is given by the formula; 

jlix) = lirn ( '''+-^^'''-- S(x -xt)) = {±1 Six ~ xt)) . (2.11) 

that is proven in Appendix A. For it to hold, it is essential to use the symmetric derivative over time of 
Xt because the left and right time derivatives lead to different results, with the difference coming from 
the white noise contribution to xt [22] . 

The probability current jl{x) may be written in the form borrowed from hydrodynamics as pt{x) vl{x) 
where 

vl{x) = pt{x)-'jl(x) = = ^*(^) - rft'(^)a.lnpt(a;) (2.12) 

[d{x- Xt)) 

has the interpretation of the time t mean velocity of the process conditioned to be at point x (once 
again, the velocity should be defined by the symmetric time derivative) . Accordingly, the quantity vt {x) 
is called the mean local velocity. Geometrically, vt is a time dependent vector field on X, as we 
show in Appendix B. The continuity equation p.9|l takes now a hydrodynamical form of the advection 
equation 

dtPt + V • [ptvt) = (2.13) 

for the density pt{x) transported by the velocity field Vt{x). 

The vanishing of the probability current jt{x) for densities pt, or of the related mean local velocity 
Vtix), is usually taken as the definition of the detailed balance for the process xt- It assures that the 
instantaneous PDF of xt is time-independent; pt = p. Assuming the detailed balance and introducing 
the Hamiltonian H{x) = — /3~^ Inp(a::) -I- const., where (5~^ is the temperature in the energy units, the 
SDE (|2.1|l may be rewritten as the equilibrium-type Langevin equation 

x\ = -l3d','{xt){d,H){xt) + rl{x) + Q{x), (2.14) 

with the notations of (|2.4p . Conversely, a dynamics governed by equation (I2.14p satisfies the detailed 
balance relative to the Gibbs density Z~^e~^ where Z (the partition function) is the normalization 
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factor. Thus the equihbrium form (I2.14|) of the dynamics is equivalent to the vanishing of the mean local 
velocity, the property independent of the choice of coordinate system. The presence of the correction 
rt in Eq. (|2.14|l assures that the drift term transforms as a vector field under a change of coordinates if 
^-I3H transforms as a density, see Eq. (IB.10|I in Appendix B. 

The above general considerations carry over, at least on an informal level, to diffusion processes in 
infinite-dimensional spaces described by stochastic partial differential equations. Nevertheless, as ex- 
plained in Sect.[6l in few important examples of infinite-dimensional non-equilibrium diffusions there are 
obstructions to the reahzation of the part of our program that we discuss in the next section. 



3 Lagrangian picture of diffusions 
3.1 Lagrangian frame of mean local velocity 

Recall that in hydrodynamics the motion of fluid particles in the Eulerian velocity field Vt (x) is described 
by the ordinary differential equation 

X = vt{x) (3.1) 

that generates the flow x t-^ ^t{x) assigning to the initial condition x of the fluid particle at time to 
its position at time t. One has: 

dt$t{x) = vtiM^)) and $to{x) = x. (3.2) 

We assume below that "l>t is well defined for all times, see, however, Sect.|6] The passage to the Lagrangian 
frame of the velocity field Vt is realized by the family of inverse transformations x i— > $^^{x) retracing 
back the flow. We have assumed that the Lagrangian and the Eulerian frames coincide at time to. 

Let us apply the above hydrodynamical idea to the diffusion process xt, describing it in the Lagrangian 
frame of the mean local velocity Vt{x). In this frame, the process Xt becomes 

it = 9^\xt). (3.3) 

In words, Xt is the point that the particle of the hypothetical fluid moving with the mean local velocity 
occupied at time to if at time t it is at xt. We shall show that the Lagrangian-frame stochastic process 
Xt is again a diffusion by finding the SDE that it obeys. 



3.2 Instantaneous densities in the Lagrangian picture 

Let us start by addressing the question what are the instantaneous PDF's of the Lagrangian-frame process 
Xt. These are defined as 

ptii) = {Six -it)) = {S{S:-<i>T'{xt))). (3.4) 

Changing variables inside the delta-function on the right hand side, we may rewrite the above relation as 
the identity 

ptix) = (pt{x) {S{$tix) - Xt)) = (y3f(i) pf($f(i)) , (3.5) 
where ipt is the Jacobian of the transformation "I>t : 

MS:) = det((9,$t)'(i)) = [det{d,<i>ty){<i>t{x))y\ (3.6) 

On the other hand, it is well known (and easy to check) that the solution of the Cauchy problem for the 
advection equation p.l3|l may be written in the form 

pt{x) = J S{x-^t{y))pto{y)dy = <^t(i)"Vto(^) ■ (3-7) 

for X = (a:). In words, Eq. (|3.7p states that pt{x) is equal to the density ptgix) at the initial point 
of the Lagrangian trajectory passing through x at time t, divided by the factor ifit{x) giving the volume 
contraction around that trajectory. Comparing Eqs. (|3.5|l and (|3.7p . we infer that 

Pt{x) = ptoix). (3.8) 

This shows that the instantaneous PDF's freeze in the Lagrangian frame to the time to value 
of the Eulerian-frame density. Since the process Xt itself is, in general, non-stationary, this might 
come as a surprise, although it is a direct consequence of the advection equation (|2.13|l . 
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3.3 Stochastic equation for the Lagrangian-frame process 

There are further surprises in the Lagrangian frame resulting in a simphfication of the non-equihbrium 
dynamics. Let us find the stochastic equation obeyed by the process Xt- This is a straightforward, al- 
though somewhat tedious, exercise. By the standard chain rule, that holds for the Stratonovich stochastic 
equations, 

i\ = {dt^T^Yixt) + {du^^^f{xt)x1. (3.9) 
Differentiating over time the identity "1>7^ (*l't(a;)) = x and setting x = ^tix), we infer the relation 

{dt^T^)\x) = ~{dk^T^r{x)v'^,{x) = -{dk'^^^nx)[u'i{x) - d1\x)di\npt{x)\. (3.10) 
The substitution of the last equality and of Eq. (12. ip to the identity (|3.9|l gives: 

i\ = {du^T^y{xt)[~u'i{xt) + df{xt)di\npt{xt) + u'tixt) + C,'l{xt)\ 

= {dk^T^y{xt)[rl{xt) + df{xt)di\npt{,xt) + C"(xO], (3.11) 



where the second equality follows from Eqs. \2A\ . Note the disappearance of the drift ut from the right 
hand side. Let us introduce the Lagrangian-frame white-noise vector field 

cKi) = (du^^^nx)c!^(x) (3.12) 

for X — ^t{x). It has mean zero and covariance 

(Cl(i)Cl(y)) = 25{t-s)Dnx,y) (3.13) 

with 

Di^x,y) = {d,<f^'yix)D^'ix,y){d,$^y{y) (3.14) 

for X — $t(i) and y — <l?t(j/). Observe that the covariances D'^ and Dj-* are related by the standard 
tensorial rule of transformation under the map $7^. We shall need two identities that may be obtained 
from the change-of- variables relations (|B.8P and (|B.9P of Appendix B if we set ^ = there. They 
are: 

f':ix) = d^,D';^{x,y)\y=, = {d,<f^y{x)[rUx) + {d,$t)^i£)d'!\x)idHdi^^yix)] (3.15) 

and 

{di^;'y{x){dj In pt)(x) = {diinpt){x) - {d,'i>tf(x)(didh'i>^'y(x). (3.16) 

Adding the first of the latter equations to the second one multiplied by {dk^t^y{x) dt'{x), we obtain 
the identity 

f';{x) + i,\x){d,\npt){x) = {du^T'y{x) [r'^ix) + df{x){ddnpt){x)\. (3.17) 

Recalling that pt = ptg for all t and defining the Lagrangian-frame Hamiltonian by the relation 

H{i:) = In ptoii) + const. (3.18) 

for an arbitrary constant, the identity (|3.17p permits to rewrite the stochastic equation (|3.1ip in the form 
of Eq. (pT4)) : 

H = -PdYix,) {d,H){xt) + f\(it) + ait) . (3.19) 

This is the main result of this section: the Lagrangian frame process xt satisfies the equilibrium 
Langevin equation with detailed balance relative to the density pto(i) = Z^^e""^^'*' that 
stays invariant in the Lagrangian frame. 

If the original process Xt is stationary with M'(a;), D^-'{x,y) and the single-time PDF p{x) time inde- 
pendent then the corresponding mean local velocity field v{x) is also time-independent. The Lagrangian- 
frame process Xt, however, is non-stationary if v does not vanish, although its single-time PDF is equal 
to p{x) and does not change in time. The stationary non-equilibrium dynamics becomes in 
the Lagrangian frame a non-stationary equilibrium one with the same invariant probability 
density. 
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Figure 1: Left: theoretical invariant PDF p{6) (blue solid line) compared to the histogram 
of 30000 time values on 1500 trajectories of the processes 6t. In the insert 
the same figure for 6t undistinguishable with bare eye from the one for xt 

Right: the same figures for the process xt obtained by the change of variables 

dm 




Figure 2: Left: mean local velocity (blue solid line, everywhere positive) as compared 
to the deterministic velocity equal to the drift term in Eq. (|3.21|l (black 
dotted line changing sign, with a repulsive and an attractive fixed points 
well visible in the blowups) 

Right: Lagrangian trajectory ^^(t) of the mean local velocity with 6l{0) ~ 

In the case where the original process has time-independent instantaneous PDF's with vanishing 
probability current, the Eulerian and the Lagrangian frame processes coincide. However, for a non- 
equilibrium Langevin dynamics 

xl = -pd^{xt)idjHt)ixt) + rlix,) + Fl{x) + Cti^) (3.20) 

with a time-dependent Hamiltonian Ht or/and an additional non-conservative force Ft that generate 
non-trivial probability current, the passage to the Lagrangian frame of mean local velocity vt recasts the 
dynamics into the equilibrium form (|3.19p with a time-independent Hamiltonian and no non-conservative 
force. The same is true for the process satisfying the equilibrium Langevin equation p.l4p but with 
non-Gibbsian instantaneous densities (relaxing to equilibrium or not). 

3.4 Examples 

3.4.1 Colloidal particle on a circle 

The simplest example of a non-equilibrium Langevin dynamics is provided by the overdamped motion of 
a particle on a circle with its angular position satisfying the stochastic equation 

dt = -{deH){6t) + F + Ct (3.21) 
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with a periodic potential H{9) = H{9 + 2tt), a constant (non- conservative) force F, and a white noise 
("t with covariance {(^t C^s) ~ 2D5{t — s). Eq. p.2ip has a stationary solution with the invariant PDF p 
given by the formula: 

e 2n 

Pie) = Z-ie-*(«W-^'')(|e*(^(^)-^^)d^ + e^ | e , (3.22) 

X 

where Z is the normalization factor. The current corresponding to this density is constant: 

j = [-{deH){9) + F-Ddg]p{9) = DZ-^{e^-l). (3.23) 

The one-dimensional dynamics becomes simpler in the variable 

e 

X = J p{^)d&. (3.24) 



taken modulo 1. Note that ^ — p{6) = jv{9)^^ so that j^^x is the time that the Lagrangian trajectory 
OL{t) of the mean local velocity starting at 6 — takes to get to 6. In the variable x, the invariant 
density p(x) = 1 and Eq. (|3.2H) takes the form 

xt = i + r(xt) + C,t(xt), (3.25) 

where C,t{x) — p{9)(,t for 9 — 9L{j~^x) and 

r{x) = Dp{9)d,p{9) = D{dep){9) = [- {d0H){9) + F]p{9) - j . (3.26) 

In the variable x, the mean local velocity v{x) = j. The corresponding Lagrangian-frame process 
Xt = Xt — j(t — to) and it satisfies the equilibrium-type Langevin equation 

it = ft{xt) + 6 (St) (3.27) 
with ft(x) = r[x + j{t — to)) and (t(x) — C^t{x + j{t — to)) and a constant Hamiltonian. 

Fig. 1 and Fig. 2 represent the invariant density and the mean local velocity with its Lagrangian 
trajectory, both for the process 9t satisfying Eq. (13.21(1 with H{9) = 0.87s~^ x sin(6'), F = 0.85s~^ and 
D = 0.036s-\ Such process models the dynamics of a colloidal particle kept by an optical tweezer on a 
nearly circular orbit in the experiment described in [11) . 

3.4.2 Linear stochastic equations 

A general class of explicitly soluble examples of non-equilibrium dynamics, with multiple applications, is 
provided by stationary linear SDEs in d dimensions of the form: 

it = Mxt + Ct , (3.28) 

where M is a matrix whose eigenvalues have negative real part and where 

{Ct Ci) = 2D'-'3it~s), (3.29) 

with a positive matrix D = {D^-'). Here, the invariant density has the Gaussian form [5] 

pix) = Z-e-''-(^) ^ (3.30) 

H{x) = ^x-C'^ for C = 2 j e'^' D e'^'^ dt . (3.31) 

The time-integral in the formula for the covariance C converges due to the assumption on the eigenvalues 
of M. The mean local velocity corresponding to p{x) is 

v{x) = {M + DC'^)x (3.32) 

so that it depends linearly on x. The Lagrangian-frame process 

^ ^_(A/+DC-i)(t-to)^^ (3 33) 

satisfies the time-dependent equilibrium-type linear Langevin equation 

it = -l3DtVH{x) + Ct (3.34) 
where the white noise Ct = e"'-^'^^^'^ ){*-*o)(*-^ y^g^ ^j-^g covariance 

(CtC) = 25it^s)Dl' with ^ o-(^'+i=c-^K^~to)jj^-(M+DC-^)^(t-t„)_ (335) 
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3.4.3 Sheared suspensions 

Stochastic equations of the type (|2.1|l may be used to model the dynamics of suspensions of colloidal 
particles [TJ or of a polymer, undergoing an overdamped motion driven by conservative forces and 
opposed by friction, see |29) for a recent discussion. An example is provided by the set of equations for 
the three-dimensional positions of A'^ particles: 

-yfi = -9r°-ff(r) - 7u"(ri) + dt , (3.36) 

where 7 is the friction coefficient, r = (ri)^i, H(r) is the potential energy and u{t, r) is the velocity 
field of the solvent. are the components of the white noise with the covariance 

{(.ttils) = 2 7/3-^5"" 5,, <5(t-s). (3.37) 
For a diluted colloidal suspension, assuming only 2-body isotropic interactions, one may take 

H{r) = Y.^{nj) + (3.38) 

for Tij = \ri — rj\ and for the polymer modeled as a chain of beads with nearest neighbor interaction 
(Rouse model [27]'). 

H{r) = + E '^»(^») • (3.39) 

i<N i 

If the solvent is at rest, and the external potential Uo is confining then the detailed balance holds for the 
normalized Gibbs density po(ll) = ^"^e^''^'-' which is left invariant under evolution. If, however, the 
solvent undergoes a shear flow with Ut{r) — f{r ■ ei) 62, where are the vectors of the canonical basis 
of R'', or a vortical motion with ut{r) = g{\r A 63]) r A 63, then the detailed balance is broken and the 
mean local velocity becomes equal in the stationary state to 

vt{r) = -7"'9,aff(r) - ut{r,) - (7/3)"' In p(r) . (3.40) 

where p(r) is the non-Gibbsian invariant density. 

In general, the form of p(r) is difficult to access in realistic situations. One of exceptions is the 
idealized case leading to the linear stochastic equations describing the simplest realization of the Rouse 
model of a polymer suspension with U{r) = ^kt'^ and Uo{r) — \kr^ and with linear velocity field Ut{r). 
For the vortical velocity Ut{r) = wr A 63, the stochastic equation (|3.36p takes the form (|3.28p with 

MiJ" = -'y~^S"'''{-KAij + kSij) +u;eab3Sij, (3.41) 

where Aij — {Siij — Sij), and with the matrix 

— i I — 1 

D^" = (7/3)-'r''5.,. (3.42) 

in the noise covariance (|3.29|l . In spite of the vortical motion of the solvent, the Gibbs density po(2l) 
independent of the vorticity uj remains invariant for the symmetry reasons. Nevertheless, for co ^ 0, the 
detailed balance is broken and the mean local velocity is given by the solvent velocity 

vt{r) = ior.Aea. (3.43) 

The Lagrangian frame just rigidly rotates around the third axis with the angular velocity $1 = lj 63 and 
the Lagrangian-frame process fi^t satisfies the stochastic equation p.36|l with ut set to zero. 

Keeping the same harmonic potentials but replacing the vortical solvent motion by the shear flow with 
Ut{r) — s (r ■ ei) 62 with a constant shear rate s, we obtain the linear stochastic equation p.28p with 

Mf/ = -7-'5'"'(-kA., +H.,) + sr'5"5., (3.44) 

and the noise covariance as before. The N x N matrix —A = (—Aij) has the eigenvalues ujk = 
2[l — cos (7^)] corresponding to the normalized eigenvectors 

(^') = ((f)''"cos(^^^)) for ^ = 0,1,...,A^-1. (3.45) 

The passage to the Fourier modes fe = r^t^j (sum over j) diagonalizes matrix M into 3x3 blocs 
with the entries M"** — S'^^'fj.i + s S"'^S''^ for fie = kui + k. The invariant density p(r ) is Gaussian. Its 
covariance depends quadratically on the shearing rate s and is composed of the 3x3 blocs 

Cf = j-K' + a^(r'5"' + r'(5") + 2a|r'5"'], (3.46) 
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Figure 3: Ellipses followed in the xy-plane under the Lagrangian flow of mean local velocity 
by the Fourier modes ri starting at (1,0) for different values of the parameter 

'^e = — 7 — r\ 

4K(l-cos(^))+2fe 



for ae = = „, with the blocs of the inverse covariance: 

(C-')f = iSfiels"'' + ^ - r^^'h - ^(^'j^^' + r'^")]. (3.47) 

L l + ^ + J 

The mean local velocity has the Fourier components 

C(r)? = {M + DC-')ffl, (3.48) 

see Eq. ((332j, with 

{M + DC-^)f = i^l (r'5" - r'^^'') - ^ (r^^"'' + ^2(5") + 2(7^^25"] (3.49) 

and is incompressible. Its Lagrangian flow is linear. It factorizes for different Fourier modes and takes 
place along ellipses in the planes orthogonal to 63: 

9t(r)t = re + (5 + 6 re) cos \ 1 + (° ^e ~ S re) — 



7i 



-(rV| + rM)^ + 25"V,yi + a|] sin(|i^), (3.50) 

see Fig. 3. The ellipses are more and more elongated in the direction of the flow with increasing shearing 
rate s and decreasing Fourier mode i. The formula for the time-dependent covariance Dt of the noise 
in the Lagrangian- frame process is given in Appendix C. 

3.5 Back to Eulerian frame 

When passing to the Lagrangian frame, a part of the information about the system contained in the 
probability current or mean local velocity is lost. If we want to reconstruct the original Eulerian process 
xt, we have to supply the forgotten information. A convenient way to do that is to provide the local 
velocity transformed to the Lagrangian frame: 

v\{x) = {dk^-^)\x)v\(^x) = -~{d,<^-,^)\x) (3.51) 

for X = ^t{x). Given the vector fleld Vt(x), consider the flow of transformations x ^ '^t(x) such that 

dt^t{x) = ~vt{^t{x)), kt^{x) = X. (3.52) 

The comparison of Eqs. (|3.51|l and (|3.52|) shows that $t(a;) = ^~[^(x). This permits to reconstruct the 
original process as 

xt = ^t^(xt) (3.53) 



9 



and the original mean local velocity as 

vl{x) = {d,^T^y{x)il[x) = -{dST^)\i) (3.54) 

for X — $t{x). In the special case when the mean local velocity is time-independent (for example when 
xt is a stationary process), 

vl{x) = ~{dt^^y{x) = -{dt^^tYix) = = v'{x) (3.55) 

so that the velocity field Vt coincides with the mean local velocity of the Eulerian frame and is time- 
independent. 

As we see, the knowledge of the non-equilibrium diffusion xt is equivalent to the knowledge of the 
equilbrium diffusion it and of the (deterministic) velocity field Vf 



4 Diffusion Processes with Hamiltonian forces 
4.1 Modified probabifity current and mean local velocity 

In many applications, one deals with non-equilibrium diffusions in the presence of Hamiltonian forces. It 
is then useful to single out their contribution and to replace the SDE (|2.ip by 

xi = ui{x) + (xO {d,H)ix,) - {d,nY){xt) + Cd^t) (4.i) 

where the term flj-' djHt with IIj-' — — IIj* stands for the Hamiltonian force. Geometrically, the an- 
tisymmetic tensor field H^-' represents a (possibly time dependent) Poisson structure but we shall not 
need its property that assures the Jacobi identity of the Poisson bracket. The subtraction of /3~^9jH(-' 
on the right hand side of Eq. (|4.1[) assures that the terms involving Ht transform as a vector field under 
a change of coordinates if the Gibbs factor e"''^ transforms as a density. An example of dynamics (14. ip 
is provided by the Langevin equation 

xl = ~l3dl'{xt) {d,Ht){xt) + riixt) + Fi{xt) 

+ TV^{x,) {d,H){xt) - (a,Hn(:,0 + C{x,) , (4.2) 

compare to Eq. (|3.20p . In the presence of Hamiltonian forces, it is convenient to redefine the probability 
current as 

jl = [u\ + K'{djH,) - h'^d,]pt, (4.3) 

where ftj"' = dj-* — l3~^Yi\^ . The new expression for the current jt differs from the one prescribed by 
Eq. (|2.10|l by the addition of the term f5~^dj{jT^ pt)- The continuity equation (|2.9|l still holds since the 
added term is divergence-less so that the flux of jt through the boundary of any region V still gives the 
rate of change of the probability that xt belongs to V. For the case of the Langevin equation (14. 2|) . the 
new expression for the current reduces to 

jl = [ - /3 d,Ht + Fl - d,] pt . (4.4) 

In particular, \i Ht = H is time-independent and the additional force Ft = Q then the modified proba- 
bility current (|4.4p associated to the Gibbs density p{x) = Z~^e~^ ^^-^^ vanishes and p is preserved by 
the evolution. It is then natural to extend the notion of equilibrium dynamics to such a case. 

As before, we may introduce the velocity field by the relation 

vl = Pt'ft = ul + nl^d.Ht) ~ bYdjlnpt. (4.5) 

Since now 

= ^/Ir^^^'u^ + p-'Mx)-'d,{U^^pt){x), (4.6) 

{S(x-Xt)) 

we shall call vt{x) the subtracted mean local velocity. The continuity equation (|2.9p still takes the form 
of the advection equation (I2.13|l . 

If we realize the passage to the Lagrangian frame of the velocity vt of Eq. (14. 5p as described in 
Sect. 13. ll using the flow of Vt that we shall still denote by $t and introducing the Lagrangian-frame 
process xt = ^^^i^t), then the considerations of Sect. l3.2l go unchanged because they only use the 
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advection equation (|2.9p . not the explicit form of Vt{x). As before, we infer that the instantaneous PDF 
of the process Xt is frozen to the time to value ptg of the PDF of the Eulerian process Xt- 

On the other hand, in the derivation of the SDE for the Lagrangian-frame process in Sect. 13.31 the 
explicit form of vt{x) was used in Eg. p.lO|l . As a consequence, the SDE for ij will pick now the 
additional term 

-idk<i>^'y{xt) r'ptixr'di (nf pt) (xt) 

= ^r\dk<S>^'y{xt)[a^'{xt){dilnpt){xt) + (anf)(a;0] 

= -r\d,<f^'nx,) [n^\x,)idi<f^'y{xt){d,pt,)ix,) 

+ ii1'{xt) {dj<i>tf{it) {didH<i>yy{xt) + {diu'^'yxt)] , (4.7) 

where the second equality follows from Eq. (I3.16P and the identity pt ~ ptg ■ Introducing the Lagrangian- 
frame antisymmetric tensor field 

uYix) = {d,<fy'y{x)u':\x)idi<i>y'y{x) (4.8) 

where x — (x) and observing that 

{d,ni'){£) = [{dMy'y{x)ii'^'{x){di'S>yy{x) + {dt<s>y'y{x){dHn1'yx){di$yy{x) 

+ id,<fy'y{x)n':'{x) idHdi<s>y'y (x)] {dj<ft f{x) 
= {d,$y'yix) (9,nr)(x) + id,<fy'y{x)u':'{x) {d,di<fyy{x) (9,<E.,)"(i) , (4.9) 

we may rewrite the additional term (14. 7|) as 

-r'fins=t){d,pt,){x,) - r\d,u^){xt) = nnxt)d,H{xt) - r^d^n^yx,). (4.io) 

Altogether, the Lagrangian-frame process Xt satisfies now the equilibrium-type time-dependent SDE with 
a Hamiltonian force: 

ii = -f3dY{xt) {d,Hyxt)dt + fiYixt) {d.nyxt) - i3-\d,&yyxt) + -f cK^O , (4.11) 

Clearly, the modified probability current associated with the conserved density pt„ = Z^^e^^^ vanishes 
for the Lagrangian-frame process. 

4.2 Example of Langevin-Kramers dynamics 

The particular case of Langevin dynamics with Hamiltonian forces is provided by the 2""^ order Langevin- 
Kramers SDE 

m^iq^ = -7»j9t - ^^Vt{qt) + U{qt) + Ct.» (4.12) 

with the positive mass m = {rriij) and friction 7 — {•yij) matrices that, for simplicity, we assume 
independent of t and q, with a potential Vt{q) and a non-conservative force ft{q), and with a white 
noise with the covariance 

{^t,i^sj) = 2/?-V,j,5(t-s). (4.13) 

We keep the matrix a different from 7 to allow noises modeling environments with variable temperature 
that violate the Einstein relation a — ^. The 2""^ order equation (|4.12p may be rewritten as the I''* order 
SDE (|4.2|l in the phase space of points x — {q,p) if we set 

d = (0 ^)-i<t) ' n = , Ht ^ hp-m-^P + Vt(q), 

Ft = (0, (a-7)m-V + /t(g)),, C* = (0, it). 

The subtracted mean local velocity in the phase space has here the form 

vt = {m-^p + p-^\Ip\npt, --VVt-7m"'p-f /t-/3"'V,lnpt-/3-VVplnpt) (4.14) 

and it vanishes for the Gibbs density p{q,p) ~ Z~^e~^ ^'-'''^^ in the equilibrium case where ct = 7, the 
potential Vt is time-independent, and the non-conservative force ft is absent. 
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4.2.1 Harmonic chain 



An example of a Langevin-Kramers dynamics is provided by a Fermi-Pasta-Ulam chain f9| with ends 
coupled to a friction force and a white noise. Such chains were often used in the theoretical studies of the 
Fourier law [2]. Here q = (rf ) with i = 1, . . . , A*', a = 1, . . . , d, and 



7?/ = 70 5'"'(<5a5i, + S.nSnj) , af/ = 70 5"''((1 + r?)5a<5ij + (1 - v)3^NSNJ) , 

JV 

mf^ = mo^''^.,, Viq) = E ^^(^(.+1)) + Et^o(rO, (4.15) 



The dynamics in the bulk (i.e. for i 7^ 1,A'') is purely Hamiltonian, whereas the boundary degrees of 
freedom ro and rjv are exposed to the thermal noise at temperatures /3~^(1 ± rj), respectively, and 
to friction. The harmonic case (that does not lead to the Fourier law 25 ) with U{r) = ^r^ and 
Uo{r) — ^r'^ corresponds to the Unear stochastic equation of the type p.28p with the matrices 



Dfj - P So<5"''(o {.-l+v)s^lSlJ+{.'^-r^)SiNSNj) ■ (4-16) 
The covariance matrix of the invariant Gaussian measure has the form 

^ p-'5^'(i-^)., (4.17) 

with matrices X, Y, Z that may be calculated exactly [25] (for rj — 0, it reduces to the covariance of the 
Gibbs measure). The subtracted mean local velocity is 

v{q,p) = (Af + DC-^-/3-'nC-^) (p, (4.18) 

where the matrix on the right hand side has, up to terms quadratic in the relative temperature difference 
77, the entries 

r]mo S ( } \ I \ 1 • (4-19) 

\-[Y{-KA+k)) -({-«A+fe)zj .^-7om(, l(iiiyij+'5,]vyN3)+7o('5ii«ij-iijv*Nj)/ 

The Lagrangian flow of w is obtained by the linear action of the matrix e'*^^^*"^ ~^ ){*-*o) ^Jijgii 

is straightforward to calculate in the linear order in 77. 



5 Fluctuation-dissipation relations 

5.1 Equilibrium Fluctuation-Dissipation Theorem 

The equilibrium Fluctuation-Dissipation Theorem [231 SI [TU] relates the spontaneous dynamical fluctua- 
tions in an equilibrium state to the relaxation dynamics after a tiny perturbation out of the equilibrium. 
It holds for a wide class of equilibrium systems including the ones described by the equilibrium Langevin 
equation 

x\ = -P<r,'{xt) {djH){xt) + Tvl^xt) {d,H){x,) - {djTvinixt) + rlixt) + Cti^t) (5.1) 

of the type discussed above. We assume that the process Xt has the time- independent Gibbs instantaneous 
PDF p{x) = Z"^e"''^'^' and denote by ( - ) the dynamical expectation. The FDT asserts that [20] 

d.{0'{x^)0\xt)) = r'^\,^,{0'\x,))^ (5.2) 

for s < t, where 0'^{x) are functions (well behaved at inflnity), that we shall call (single-time) observ- 
ables, and where on the right hand side the expectation ( — )^ involves the process obtained by replacing 
the Hamiltonian H{x) in the original dynamics (|5.ip by its time-dependent perturbation H{x) — ht 0^{x) 
within some time interval. The left hand side is the time derivative of the 2-time correlation function 
in the dynamics determined by Eq. (15. ip and the right hand side is the response of the single-time cor- 
relation function to a small dynamical perturbation of the Hamiltonian of the system. The temperature 
/3~^ appears as the coefflcient relating the two functions. For the sake of completeness, we give a proof 
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of the FDT (|5.2|) in Appendix D. It is often more convenient to consider the time-integrated version of 
the FDT: 

{0\x,)0\x,)) - {0\x.)0\x,)) = /?-^j,^^^^(O^M),„„, (5.3) 

where ( — )^ ^ corresponds to the expectation where the original Hamiltonian H{x) is replaced starting 
at time s < t by its time-independent perturbation H{x) — hQO^{x). 

5.2 Modified Fluctuation Dissipation Theorem 

We may immediately apply the FDT to the Lagrangian-frame process xt obtained from the process xt 
satisfying the Langevin equation (|4.2p . Indeed, as was shown in Sect . 14.11 the process xt = 't^^{xt), 
where "I>t is the flow of the subtracted mean local velocity (|4.5|l . satisfies the equilibrium stochastic 
equation (|4.11|l and has the time-independent instantaneous PDF pto(i) = Z~^e~^ We infer that 

for observables 0'^{x), 

d.{d\xs)d\x,)) = /?-i^|-^,,(6^(iO),; (5.4) 

where ( — )^ involves the process obtained by replacing the Hamiltonian H{x) in the Lagrangian-frame 
dynamics (|4.11|) by its time-dependent perturbation H{x) — htO^{x) during a time interval. Observe 
that this perturbation corresponds to the replacement of the Hamiltonian Ht{x) in the original equation 
(|4.2[) for Xt by Ht{x) — htO^ {^^^ (x)). Indeed, the latter replacement adds the term 



(3htb]^{xt)d,,d'{<l>t\x))\^=^, (5.5) 
on the right hand side of Eq. (|4.2|l and, in virtue of Eq. (|3.9p , results in the additional term 

{d,<ft'y{xt) [l3ht 6r(^t)a,,U..,6^(<&7'(^))] 

= I3htbl'{xt){d,d'){xt) (5.6) 

with b]-' — — f3~^lV^ in Eq. (|4.1ip for xt — $7^(^t) (with the same transformations <l>t as in the 
unperturbed process). Upon deflning the Eulerian-frame time-dependent observables 

Onx) = d^{^T\x)), (5.7) 

the Lagrangian-frame FDT (|5.4|l may be rewritten as the identity 



ds{Ol{x^)0^t{xt)) = /3-i-|-|-^^^(o?(a;0)^. (5.8) 

Note that the time-dependent observables Ot (x) are constant along the Lagrangian trajectories of the 
velocity (|4.5|l : 0"("l>i(a;)) = 0°'{x). In other words, they obey the scalar advection equation 

dtOt + vt ■ VO" = (5.9) 

and are frozen in the Lagrangian frame of the subtracted mean local velocity Vt- Since the values of the 
time-dependent observable may be chosen arbitrarily at time s and that of O'^ at time t, the only 
trace of time dependence of the observables 0° in the identity (|5.8p for flxed pair of times s < t enters 
through the time derivative ds on the left hand side that differentiates also the explicit time-dependence 
of determined by Eq. (|5.9p . We may then rewrite Eq. (15.8(1 using observables frozen in the Eulerian 
frame as the Modified Fluctuation-Dissipation Theorem, 

d.{0'{x,)0\xt)) ~ {{v,-VO'){xs)0\xt))^ r'sh\h=o(Oti^t))^, (5.10) 

where the expectation ( — ) ^ on the right hand side refers now to the process obtained by replacing the 
Hamiltonian Ht in Eq. (|4.2p by Ht{x) — ht 0^{x). In the time-integrated form, Eq. (|5.10p becomes 

t 

{0\xt)0^{xt)) - {0\x,)0^xt)) - j {{v„-VO^){x^)0^{xt))du 

s 

= /3"'lrLo=o(o'(-*)).o..' (5.11) 
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Figure 4: Left: the bottom coinciding curves: LHS (continuous black) and RHS (daslied red) 
of the integrated MFDT (fSTTj) for 0"(6l) = sin(6l), the upper (dot-dashed 
green) curve: the first two terms on its LHS, the middle (dashed blue) 
curve: the corrective integral term 

Right: RHS (dashed red curve) and LHS (continuous black curve) of the integrated- 
in-time Lagrangian- frame FDT (|5.4|l with 0'^{6) = sm{9) 



with a corrective integral term with respect to Eq. (|5.3p . An experimental check of the time-integrated 
MFDT for a colloidal particle has been described in [TT] . Fig. 4 shows the numerical check of this 
relation and of its Lagrangian-frame counterpart for the stationary process solving the SDE (I3.21|l for 
0"(6() = sin{e). 

The MFDT was proven directly in ^ in the stationary setup and shown to be equivalent to identity 
(|5.8|) similar to the equilibrium FDT (15. 2|) but for observables frozen in the Lagrangian frame of mean 
local velocity. In the present paper, we unravel the deeper reason for that equivalence, namely the fact 
that the non-equilibrium diffusion process observed in the Lagrangian frame of the (subtracted) mean 
local velocity evolves according to an equilibrium dynamics with a time-independent instantaneous PDF. 

5.3 Links with fluctuation relations 

Ref. [6]) also discussed fluctuation relation extending the MFDT to non- stationary situations. It was 
shown there that the Hatano-Sasa version [13] of the Jarzynski equality [151 116j reduces close to station- 
arity to the MFDT for special observables and that one need Croooks' extention of the Jarzynski- 
Hatano-Sasa equality to extract at stationarity the MFDT for general observables. The results of the 
present paper permit to propose yet another extension of the MFDT. 

For the process xt evolving accordingly to the Langevin equation (|2.14|) but with a time-dependent 
Hamiltonian Ht{x), the Jarzynski equality reads 

(e-'^^'o.) = |l, (5.12) 

where *° 

t 

Wto,t ^ J{dsH){xs)ds and Zt = je-^"'^''Ux, (5.13) 
to 

provided that the PDF of xt^ is Z^^^e~^^*f> . Applied to case with the Hamiltonian Ht{x) = H{x) — 
ht 0''{x) and expanded to the second order in functions /i", Eq. (|5.12p reduces to the FDT (|5.2p . 

a = l,2 

The proof goes as in [5] where it was written for a less general case. 

The above observations apply to the case of the Lagrangian-frame dynamics. For the process xt 
satisfying the SDE (|4.11|l but with the Hamiltonian H replaced by H{x) + H[{x) with H[ = for 
t < to, we have for t > to the Lagrangian-frame version of the Jarzynski equality: 

t 

^e-/3m„,^ ^ ^ ^^^^^ ^ j^Q^H[){is)ds (5.14) 

to 

and Zt = j e-'^'^"+"t'>'-^'>dx, provided that the PDF of Xtg is Z-^e'l^" = ptg. The process Xt 
such that Xt = ^t^i^t), with ^t standing for the Lagrangian flow of the mean local velocity vt{x) of 
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the unperturbed process Xt, satisfies the SDE (|4.2|l with the original Hamiltonian Ht{x) replaced by 
Ht{x) + H[{^'[^ (x)). This follows by the same argument as around Eqs. (|5.5|l and H5.6|) . In terms of the 
perturbed process xt, 

t 

Wt„,t = j{dsH's){^7\xs))ds and ^ = j "'^'^''^ pt„{x) dx . (5.15) 
to 

For H'{x) ^ - Yl htd"-{x), one obtains the Lagrangian-frame FDT equivalent to the MFDT 

a=l,2 

1)5. lOp by expanding the identity (|5.14|l to the 2"^* order in h^. Not very surprisingly, there exist different 
fluctuation relations that may be viewed as an extension of the MFDT to more general situations. 



6 Non-equilibrium diffusions without Lagrangian picture 

In the preceding sections, we have discussed difi'usion processes in a finite dimensional phase space. The 
basic assumption underlying the discussion of the Lagrangian-frame picture of diffusions was the existence 
of the Lagrangian flow x i— > ^t{x) of the mean local velocity satisfying Eqs. (|3.2p . This is guaranteed if 
the velocity Vt{x) is smooth and the (phase-)space X is compact, like in the circle example, but may be 
not assured if X is unbounded in which case the Lagrangian trajectories of vt may blow up in finite time. 
The idea of the decoupling of probability flux by the passage to the Lagrangian frame of the mean local 
velocity can, in principal, be applied to infinite-dimensional diffusive processes. It appears, however, that 
a number of known examples of diffusive processes described by stocheistic PDEs do not allow a global 
flow of mean local velocity and, hence, do not admit a Lagrangian-frame equilibrium-like description. Let 
us illustrate this phenomenon in specific cases. 



6.1 One-dimensional Kardar-Parisi-Zhang equation 

The KPZ stochastic PDE [I7j describes the fiuctuations of a d-dimensional interface with the height 
function htix). It has the form 

dtht{x) = vV^ht{x) + l\{SJht{x)f + Ct(a;) (6.1) 

where C,t{x) is the white noise with the covariance 

{U^)Uy)) = 2D S{t ^s)5{x^y). (6.2) 

The adjoint generator of the process ht in the (infinite-dimensional) space of the height functions h has 
the form 

Lt . /^[- .V^Hix) - |(V.(.))^ + ^^l'^- (6-3) 

A straightforward (although somewhat formal) calculation [12) shows that in one space-dimension with 
periodic boundary conditions (where Vft = dxh), the Gaussian density in the space of height functions 

p[h] = ^-ie-2T7/(^'^("»'''" (6.4) 

is annihilated by (for all values of A) and thus stays invariant. The corresponding mean local velocity 
given by Eq. (|2.12p has the form 

v[h\{x) = \\{Vht{x)f (6.5) 

and the Lagrangian trajectories of v[h] should be solutions of the equation 

dtht = \\(Vht{x)f (6.6) 
that becomes for Utix) = —W/ht(x) the inviscid Burgers equation [3] 

dtut{x) + ut{x)Vut(x) = (6.7) 
with the solutions satisfying the relation 

Ut{x + {t — ti))uta(x)) = Uta{x) (6.8) 

and developing discontinuities (shocks) for the first time ts > to such that ts = to -\ 3:2-3^1 — ^ f^j. 

a pair of points {xi,X2)- The corresponding height function ht{x) looses at t = ts the differentiability 
and, although weak solutions of the inviscid Burgers equation exist beyond the time ts , there is no unique 
global invertible Lagrangian flow of the mean local velocity v[h] and no global Lagrangian-frame picture 
of the KPZ evolution. 
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6.2 Diffusive hydro dynamical limits 

Similar problems obstruct the existence of the Lagrangian picture in the effective equations describing 
the large-deviations regime of fluctuations around diffusive hydrodynamical limits of some lattice particle 
systems. The evolution of the particles consists of random jumps to nearby sites. On the scales of the 
order of the size of the system L, and for times of the order L^, such stochastic evolution gives rise to 
an effective diffusion in the space of macroscopic densities nt{x) [301 118) . The dynamics of the densities 
is given by the continuity equation dtut + V ■ jt =0 for 

jlix) = ~^D^'{n,ix))d,ntix) + Ct{x\nt) (6.9) 

where ^t{x\n) is the density-dependent white noise in time and space with the covariance 

{C{^\n) Ciiyln)) = e5{t~s)5{x-y)x''{n{x)), (6.10) 

where oc L~'* is the total number of microscopic particles assumed to be large. In particular, in the 
limit where e = 0, the density nt{x) satisfies the deterministic hydrodynamical-limit diffusion equation 

dtnt = '^d,{D''{nt{x))djnt{x)) . (6.11) 

One considers such systems with periodic boundary conditions or with Dirichlet ones where one fixes 
the boundary values of the density nt{x) on the boundary of a finite domain A C R'*. The first case 
corresponds to an equilibrium evolution whereas the second one to a non-equilibrium boundary-driven 
one. The adjoint generator of the process rit has the form 



up to the terms of higher orders in e. To the leading order, the stationary PDF in the space of density 
functions takes the semi-classical form 

p[n] = e-^s["l (6.13) 
with the functional S[n] satisfying the Hamilton- Jacobi equation 

and a certain stability condition 1 . According to Eq. (I2.12p . the mean local velocity in the space of 
densities has the form 

vln]{x) = '-[d,{D^'in,ix))dMx)) ' 5^1))] (^.15) 

up to terms that vanish at e = 0. The functional S[n\ is explicitly known in few boundary driven 
non-equilibrium situations for which one may study the existence of the Lagrangian trajectories of v[h]. 

6.2.1 Zero range processes 

Here, D''^ {n{x)) = tp' (n{x)) 5^-' and x'"' ('^•(s^)) — ip{n{x)) 5'^ for an increasing function ifi > of n > 
related explicitly to the jump rates of the zero-range particle dynamics [18) . The hydrodynamical-limit 
equation ()6.11[l reduces to the form 

dtnt{x) = |vV(w(a:))- (6.16) 



and the functional S[n] satisfies the relation [T] 

5S 



Sn{x) 



lnip{n{x)) - InA(x) , (6.17) 



where X{x) = Lp{n{x)), with n{x) providing the stationary solution of Eq. ()6.16p so that X{x) is a 
harmonic function on the domain A with prescribed boundary values. In virtue of Eq. (|6.17p . 

a,(x'^(n(x))9,^) = V-^{n{x))V[ln^{n{x))-ln\{x)] 

= \/^^{n{x)) - V- (<^(n(a;)) VlnA(a;)) . (6.18) 
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One infers that in this case 

v[n]{x) = ■ {(p{n{x))V\n\{x)) . (6.19) 
The equation for the Lagrangian trajectories of v[n] has the form 

d^mix) = ^^'{n{x)) (Vn(x)) . - i^(n(^)) ^^^^ (6.20) 

which is a quasi-hnear f'-order PDE whose solutions may be composed from characteristic curves. The 
existence of global solutions will again be obstructed by caustics, i.e. by crossings of the projection of the 
characteristics to the space. That this phenomenon takes really place may be easily seen in one dimension 
where )^{x) is a linear function. 

6.2.2 Symmetric simple exclusion process (SSEP) 

Here D'-' = S^-' and x'^('^) = '^(l ~ The functional S[n] is explicitly known in one space-dimension 
[S]. It satisfies the identity [I] 

= In- — — <fi[x\n) (6.21) 



Sn{x) 1 — n{x) 

where ip{x\n) is the solution of the ordinary differential equation 



+ = nix) (6.22) 



(y<p)^{x) i + c-^t^) 

with prescribed boundary values. The mean local velocity has the form 

v[n]{x) = ^V{n{x){l-n{x))Vip{x\n)) . (6.23) 

We do not know if there are obstructions to the existence of the corresponding Lagrangian flow. 

7 Conclusions 

We have shown that non-equilibrium Markov diffusions become equilibrium ones when viewed in the 
Lagrangian frame of their mean local velocity. More exactly, the diff'usion process transformed to that 
frame, although in general non-stationary, satisfies the detailed balance and has instantaneous probability 
density that does not change in time and is equal to the Eulerian invariant density if the original process 
is stationary. The passage to the Lagrangian frame decouples the non-zero probability current from the 
non- equilibrium process. The equilibrium nature of the Langevin-frame process explains on a deeper 
level the equilibrium-like fiuctuation-dissipation relations observed in the Lagrangian-frame of mean local 
velocity in [281 [6] . Our analysis indicates that the equilibrium and non-equilibrium diffusions are closer 
than usually perceived and the entire difference between them may be encoded in the probability current 
that does not vanish in the non-equilibrium case. This seems to be an interesting observation on the 
fundamental level. In practice, although the passage to the Lagrangian frame may be realized numerically 
in simulations of small systems, its experimental realization is far from obvious and its use in the analysis 
of stationary non-equilibrium dynamics may be hampered by the absence of knowledge of the invariant 
measure that enters the expression for the mean local velocity. As we have also seen, our arguments 
apply only to diffusive systems with the global fiow of mean local velocity. Such global fiow is absent in 
important examples of non-equilibrium diffusions described by stochastic partial differential equations. 
It remains to be seen to what extent a similar analysis may be carried through for other models of 
non-equilibrium dynamics. 



Appendix A 

We check here the formula (|2.1ip for the probability current (I2.10p . First note that for a similar average 
as in Eq. (|2.11|l but with the right time derivative, 

{xl.5(x~xt)) = lim ( i^iii^ ^ S(x — xt)) 

= lim ^pt{x)(^ P(t,x;t + e,y)y' dy - x'-'j = pt{x) {Ltx^) 
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= [{^,{x) + (d,cr,'){x)\ pt(x) . (A.l) 
On the other hand, for the left time derivative, 

( x\_ 5{x — Xt) ) = hm ^ — 5{x — xt) ^ 

= lim ^ {pt{x)x^ ~ j pt-e(y)v Pit - e,y;t,x)dy^ = {Llpt){x)x' - Ll{pt{x)x') 

= [ul{x) - {d,d^){x) - 2dn^)d,] pt{x) , (A.2) 

where the second equaUty combined the derivatives over e of pt-^ and of P{t — e, y; t, x). The addition 
of the relations (|A.2p and (|A.1|) gives the identity pTTTj) . 



Appendix B 

Let us check that under the change of variables x i — > x' — 'i'{x), the mean local velocity (|2.12|l transforms 
as a vector field. In new variables, the process x'f = 9{xt) satisfies the Stratonovich stochastic equation 

x'^ =u';{x') + (B.l) 

with 

u':{x') = {dk<l'nx)u'[{x), Cr(x') = (9.*)'(x)Cf(a;) (B.2) 
for x' = 'I' (a;). The covariance of the white noise Cti^') is 

{(:n^')aiy')) = 2S{t-s)D';'ix',y') (B.3) 

for 

D':^{x',y') = (9,*)'(:r)A''(^,y)(ft*y(y) (B.4) 
and x' — 'I'(x), y' = '^{y). The instantaneous PDF of the process is 

P'.ix') = p.ix){^y\ (B.5) 

where ^^^j^^^^^ stands for the Jacobian of the change of variables. In the new variables, the mean local 
velocity ((?T2|) is 

v"{x') = Ut{x')-dt^x'){dj\np't){x'), (B.6) 



where 

d'^-' [x) — D'^-' [x ,x) and Uf{x') = Ut{x) — r'f{x) . (B.7) 
The deterministic correction 

r';{x') = dy„D';\x',yX,^^, = [d,<i'-^)\y')d^A{d^<i'y{x)Df{x,y){di^y{y)\\ 

L J I y — x 

= {dk-^nx)diDf{x,y)\y=^ + {d,-^-^f{x') (9fe*)'(x)df (x) {dudi-^y{x) 

= (9fe*)'(x) [rf (x) + df{x) {d,^-^f{x') {dudi-^Y {x)\. (B.8) 

On the other hand, using the standard formula for the derivative of the logarithm of a determinant, we 
obtain 

{di<fY{x){d,\np[){x') = [^,^,y{:c){^,<i~')\x')^u[\np^{x)-\u^^^^] 

= {di\npt){x) - {d,-^-^)\x'){didH-^y{x). (B.9) 

Hence 

rf(x') + d':^{x'){d,\up[){x') = {du^y{x)[r';{x) + df{x){d,<i'-^f{x'){dudi<i'y{x)\ 



+ 



{dkm\x)dt\x)[{di\npt){x) - {d,<i!-^f{x'){didh^y{x)\ 

= [du-^Yix) [r,'-(x) + d1\x) (di Inpt^x)] . (B.IO) 
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Finally, using also the 1''* of Eqs. (|B.2p . we obtain the identity 

vt{x') = u't (x) ~ r't{x') - d't^{x') {d J In p't){x') 

= {d,^y{x)[n'l{x) ~ rUx) ^ d':'{x){dilnpt){x)] = {dk<S>y (x) v^^) , (B.ll) 
which was to be shown. 

Appendix C 

We give here the explicit formula for the time-dependent noise covariance Dt of the Lagrangian-frame 
process corresponding to the harmonic Rouse polymer in linear shearing flow considered in Sect. 13. 431 
keeping the notations of that section. Dt is composed of 3 x 3 diagonal Fourier blocs 

+Aoj sin' ( ^L^)] + iS'^'S"^ - sin ( ^) - 2a, sin^ (^^)] } (C.l) 

that are positive matrices with constant determinant equal to (7/3) 



Appendix D 

We give here a proof of the FDT (|5.2|) around the non-stationary equilibrium dynamics described by the 
Langevin equation (|5.1|l . On the one hand, the two-time dynamical correlation function is 

{0\xs)0'{xt)) ^ J p{x)0\x)P{s,x-t,y)0\y)dxdy (D.l) 

where p{x) = Z~^e~'^ is the Gibbs instantaneous PDF of the process Xt satisfying the SDE (|5.1|l and 
P{s, x;t,y) are the transition PDF's. Using the first of the Kolmogorov equations (|2.8|l and integrating 
by parts, we infer that 

d,{0\x,)0\xt)) = -J {LlpO^)(x)P{s,x-t,y)0\y)dxdy, (D.2) 

where 

= [-pdl'idjH) + nl'idjH) - p-'{djni:^)]d, + d.dl'dj (D.3) 

and LIp — 0. Let 

= + [pdl' {d,0^) - nl' {d,0^)]d, , (D.4) 

be the generators of the process obtained by the replacement i/ — /is O^. Clearly, (Lj)t(pe'"'='^ ) = 

0. Expanded to the first order in hs, the latter equality implies that 

/3Li(pOi) = -g|-|^^JL^)tp. (D.5) 

As a consequence, 

ds{0\xs)0'(xt)) = P'' 1 {^\^^^{L'!yp){x)P{s,x-Xy)0\y)dxdy 

= I p{x){^\^^^^L':)(x)P{s,x;t,y)0'iy)dxdy. (D.6) 

The right hand side is equal to ( 0^(a;t) )^ so that the identity (|5.2p follows. 
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